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1. INTRODUCTION 

The use and supply of fossil fuels have been declining over the last ten years due to their adverse 
effects on the environment, which in turn raises the demand for renewable energy sources [1]. Support for the 
use of renewable energy sources has been enhanced, not by improving power efficiency, but by maintaining 
the climate target criteria [2] and thus increasing the dependency of the renewable energy source as a means 
of mitigating the traditional way of generating energy on the environment. Electricity is not only an essential 
factor for human society's growth and advancement but also plays a crucial role in economic and product 
production in the industrial sector [3]. Because of that, power should be available to any corner of the planet. 
One way of generating electricity is to use wind as a source of renewable energy and use it through a wind 
turbine that can convert the kinetic energy in the wind into meaningful electricity [4]. Another way is to use 
the photovoltaic cell, where the production of power is carried out by converting the solar radiation through 
the photoelectrical effect process. Despite the more outstanding performance and cost-effectiveness of 
photovoltaic panel technology, the use of wind power has increased in recent years. In 2017, for example, 
34% of wind turbine growth was observed in Europe in comparison to 2016 [5]. Wind turbine technology has 
become a popular field of research because of increasing demand. The usage of wind turbines is not only 
fulfilling the need of the energy industry but also becomes an important area of researching the various 
academic field. Nowadays researchers are focusing on developing micro wind turbines installed in various 
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applications like in bicycle as a portable charger [6], in-home, and building as a source of electricity [7]. 

Various control techniques have been conducted to increase power efficiency [8]-[17] and monitor 

aeroelasticity structural damage [18], [19]. Also, many research projects have focused on intensifying the use 

of wind turbines in a turbulence environment and because of this, computational techniques have been 
utilized in wind power production platform that involves significant research areas such as micro-locations, 
wind simulation, forecasting, and blade optimization [20]. Few researchers have used the Blade element 
method (BEM) code in wind turbine applications to operate to its maximum power coefficients [21]. They 
found an ideal rotational speed that provides optimal power for a specified streamlining speed. In addition to 
it, an updated BEM theory has also been used to simulate forces on large-scale wind turbines to simulate the 
dynamic model. It has been found that this dynamic model was most appropriate for engineering purposes 

[22]. Due to the rise in demand for BEM its application is increased from designing a horizontal axis wind 

turbine blade to developing a mathematical model [23]. However, the BEM method failed to simulate the 

whole flow field. Due to this drawback of the BEM method, several investigators have employed the 
computational fluid dynamics (CFD) technique in modeling the National Renewable Energy Laboratory 

(NREL) turbines by solving Navier stocks equations [24]-[26]. The wind turbine is composed of various 

parts like blades, tower, gearbox, generator, controller and many more things and every part have their 

importance but among all other parts, blade plays an important role in determining the efficiency of the wind 
turbine. So, for an efficient wind turbine, a study on the blade is necessary to get the optimum results, and for 
this, the idea of twist angle and pitch angle is necessary. Researchers used the un-twisted blade to get the best 
angle of attack by conducting several simulations for various wind speeds. The test parameters included five- 
pitch angles and four wind speeds. The results obtained from the CFD Analysis were compared with the 
experimental results from the NREL [27]. Moreover, few authors also have examined NREL Phase VI small- 
sized wind turbines for comprehending their aerodynamic behavior using the CFD technique. To execute the 
analysis, they considered five wind speeds and a constant pitch angle and for modeling turbulence, the Shear 
stress transport (SST) k-w model was considered. They found the stall at the blade root at 7 m/s [28]. Some 
of the researchers also used code called CFD Ship-Iowa 4.5, which is dynamic and uses the property of 
incompressible to check the performance parameter of NREL phase VI. Using these code two objectives of 
the NREL phase VI wind turbine are completed. One is performing analysis by fixing pitch angle and 
varying wind speed and the second is performing analysis by fixing wind speed and varying pitch angle. In 
both the case, rotor rotational speed has been taken as 7.53 rad/sec. Here detached eddy simulation has been 
used as a turbulence model. The experimental results validated various performance parameters, such as 

power, thrust, and pressure variations around the airfoil [29]. 

As proper monitoring of the blade pitch mechanism is required, there have also been a few 
suggestions for minimizing wear in the wind turbine pitch drive [30]. Also, several valued logics such as 
fuzzy logic is utilized and optimization processes such as genetic algorithms, meta-heurestic, cuckoo search 
algorithm and Particle Swarm Optimization (PSO) have been used to monitor the pitch angle by 
accommodating nonlinearity, which further leads to stability in the wind turbine system by reducing charging 
effects on a wind turbine blade [31]-[40]. Hence from the previous literature, it is clear that pitch angle plays 
a vital role in determining the wind turbine performance. So, the purpose of this research is to examine the 
influence of a pitch angle on the performance factor of the horizontal axis wind turbine blade by considering 
different pitch angles and rotational speed at different wind speeds. This research can give the researcher data 
to design and optimize the blade effectively by identifying the optimal pitch angle for the corresponding 
speed. For the numerical simulation of the wind turbine, aerodynamic analysis for different pitch angle and 
wind speed was carried out using the Ansys Fluent 15 CFD simulation [41]. 

The novelty of this paper is described as follows: 

e In the previous literature, either lower-fidelity modeling or a few pitch angle values were investigated for 
its impact on HAWT efficiency. To better understand the effectiveness of pitch angle on turbine output it 
is necessary to consider more pitch angles, as a result, a complete understanding of turbine performance 
is lacking with the involvement of few pitch angles. The present research is complementary to the 
existing knowledge by using highly precise CFD measurements and the association of more pitch angles. 

e The working conditions examined include seven-pitch angles, two different rotational speeds, and 
seventeen different wind speeds. In previous literature, this kind of working environment was not used to 
provide a more comprehensive view of how the aerodynamic efficiency varies at the next pitch angles. 

This paper is organized as follows: Section 2 represents a methodology, which comprises three- 
component, 1.e., wind turbine model, domain and meshing, and physics setup. Section 3 presents results and 
discussions followed by the conclusion in section 4. 
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2. METHOD 

The method is divided into three parts, first is the selection of airfoil and blade design for the wind 
turbine. Second is the formation of computational domain and meshing, and the final is the selection of the 
algorithm for physics setup. 


2.1. Wind turbine model 

Figure 1(a) displays the wind turbine model considered in this article. The three-dimensional wind 
turbine blade is constructed with the geometrical parameters suggested in the report [42]. This blade includes 
3 type of airfoils that are mounted from root to tip shown in Figure 1(b). The blade and rotor parameters are 
defined in Table 1. For an effective design of the wind turbine rotor, it is mandatory to know the structural 
and aerodynamic requirements. From an aerodynamic standpoint, a thin airfoil must be preferred to get a 
higher lift. From a structural point of view, dense airfoils must be chosen for high rigidity, and failures due to 
bending must be eliminated. In the design process, meeting these two requirements is of utmost importance. 
To satisfy both aerodynamic and structural demands, various airfoils must be mounted at a separate section 
of the blade. A thick airfoil is used on the root of this paper and a thin airfoil on the tip section 1s mounted. 
This article takes one meter additionally to show the blade attached to a hub and compensates for the blade 
root by | meter from the rotation axis. This paper does not involve the hub. 
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Figure 1. (a) 3D wind turbine blade, (b) NREL airfoils 


Table 1. Important specifications of the wind turbine 


Important factors value Unit 
Power 1.5 MW 
Number of blades NB 3 Not Applicable 
Rotor Radius R 43.25 Metre 
Rotational velocities 1.57,2.22 rad/sec 
Pitch angle 0,5,10,15,20,25,28 Degree 
Velocities 8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24 m/s 


2.2. Domain and meshing 

The model will use a periodicity boundary condition. ANSYS builds equations that make the 
solution on the 0° and 120° planes equal. When the results are expanded, the 240° plane will have the same 
values as the 0° and 120° planes. The values on the 60° plane will be the same as the values on the 180° 
plane after the results are expanded. The width between the domain inlet and outlet is 270 m. The global 
origin is placed in the middle of the blade root when constructing the domain. The distance from the blade to 
the velocity inlet and the pressure outlet is shown in Figure 2. The upper portion of the computational domain 
is considered as top velocity and is given the same input as inlet velocity. The radius at inlet and outlet s 
taken as 120° degree and 240° degree. The computational domain and boundary conditions for the wind 
turbine model are depicted in Figure 2. This paper uses structured grids which are symbolized by tetrahedral 
elements in 3D. Thus, the mesh of the fluid domain reflects tetrahedral form, as seen in Figure 3. Also, the 
prismatic inflation layer is added for outward flow to trap the boundary layer on an entity. The sphere of 
influence was used with a radius of 30 m and an element size of 2m to establish further precision along the 
circumference. 
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Figure 2. Computational domain with boundary Figure 3. Mesh of the computational domain 
conditions 


2.3. Physics setup 

In this paper, a pressure-based steady-state CFD simulation has been performed. The turbulence 
model used for this study is k-omega SST. In solution methods, pressure-velocity coupling was implemented, 
and algorithms used to calculate gradients and pressure are the least-squares cell-based algorithm and the 
standard algorithm. Moreover, for evaluating the momentum equation second-order upwind algorithm and 
for turbulent kinetic energy, and specific dissipation rate first-order upwind algorithm were taken. The 
residuals are monitored to 1500 iteration to reach the convergence. 


3. RESULTS AND DISCUSSIONS 
3.1. Estimation of torque 

As the blade is generally designed for many pitches and wind velocities, the optimum pitch angle 
will be activated to an automated control mechanism for a particular wind velocity to generate maximum 
torque. Torque generation by the blade not only depends on the velocity of incoming airflow, but it also 
depends on the rate of change of momentum, as the frontal area changes for different pitches, which in result 
gives the modified rate of flow and leads to torque generation. Here Figure 4 (a) and Figure 4(b) shows that 
when the pitch angle 25° and 28° are considered, and the flow velocity varies from 8 m/s to 24 m/s, torque 
decreases monotonically. In contrast to other pitch angles initially it falls and then increases. The maximum 
torque obtained is corresponding to a 28° pitch angle for a wind velocity of 8 m/s. From Figure s 4(a) and 
Figure 4(b), it has been observed that as the rotational velocity increases, torque generation by a wind turbine 
rises corresponding to all pitches and wind velocities. It is to be noted when m=1.57 rad/sec and pitch angle 
0° is considered for the analysis, the graph of torque shows an increasing trend. But, when @=2.22 rad/sec 
and a pitch angle 0° is considered there is a fall in torque value when velocity changes from 8m/s to 10 m/s, 
and after that, it rises with an increase in wind speed. 
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Figure 4 Torque variation curve at different pitch angle for (a) œ = 1.57 rad/sec, (b) œ = 2.22 rad/sec 


3.2. Estimation of power 
The numerical power can be easily calculated using the following equation after the torque value is 
received from the CFD simulation. 
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where, P = power (MW) 

T = torque (N — m) 

w = angular velocity (rad/sec) 

The variation of the power curve shown in Figure 5(a) and Figure 5(b) is the same as Figure 4(a) 
and Figure 4(b). Here in Figure 5(a), pitch angle 25° and 28° shows the decreasing trend, and it happens in 
such a way that from velocity 8 m/s to 16 m/s, the curve shows flat fall, but after 16 m/s, it decreased sharply. 
But in Figure 5(b), when the blade is given pitch angle 25° and 28” the curve shows flat decrement 
consistently irrespective of wind speed. In both the above cases, despite the decreasing trend, it has been 
observed that power is coming to be more than expected for 25 and 28° pitch angles, which are more than the 
rated power output value shown in table 1, it is showing because losses like mechanical losses, rotor control 
system losses, pitch control mechanism losses, gear & shaft losses, etc. are not considered while performing 
the simulations. The variation of pitch angle 25 and pitch angle 28 is shown in Figure 5 (a) and the variation 
of pitch angle 20° shown in Figure 5(b) resembles each other. The reason behind the variation in the power 
curve is due to stall phenomena. 





. 55000 
= e gano S pich 5 
itc m 45 —e Pitch 
_.— Pitch 10 = 40000 tons 
—-— Pitch 15 x 35000 -à Pitch 10 
—— Pitch 20 i 30C? a Pitch 15 
—~— Pitch 25 z 20000 + Pitch 20 
~ Pitch 28 Q Too —— Pitch 25 
5 10 15 20 25 5 
Wind Speed (m/s) Wind Speed (m/s) 
(a) (b) 


Figure 5 Power variation curve at different pitch angle for (a) œ = 1.57 rad/sec, (b) œ = 2.22 rad/sec 


3.3. Estimation of power coefficient (Cp) 
The power coefficient can be calculated using (2) given below. 





TXW 
p= 0.5pAV3 (2) 


where Cp = Power coef ficient 

T = torque (Nm) 

w = Angular velocity 

A = Area of the rotor 

V = wind speed (m/s) 

The power coefficient is an essential non-dimensional parameter in determining the performance of 
a wind turbine. The power coefficient curve is shown in Figure 6(a) and Figure 6(b) gives an idea about 
which pitch angle is best to work for and which pitch angle can lead to random behavior. According to the 
previous studies, the maximum value of power coefficient is set to 0.59 according to Betz's law [43], so 
keeping this limit in consideration and visualizing Figure 6(a) and Figure 6(b) one can say that values which 
are more than 0.59 are practically not feasible but in CFD the reason behind of getting such a big value is that 
losses like mechanical losses, rotor control system losses, pitch control mechanism losses, gear & shaft 
losses, etc. are not considered while performing the simulations. From Figure 6(a) it has been observed that 
for @=1.57 rad/sec when the blade is given pitch angle 25° and 28° best power coefficient obtained at a 
velocity 10 m/s, 1.e., 0.56, which is closer to Betz's limit. Blade with pitch angle 20° also reflects the best 
power coefficient value as 0.52 when wind speed is considered as 9 m/s. blade with pitch angle 15 and 
velocity 8 m/s and blade with pitch angle 28° and velocity 11 m/s also shown the good results of the power 
coefficient as 0.47 and 0.41. Figure 6(b), on the other hand, shows some different variations of power 
coefficient due to different angular velocity, 1.e., 2.22 rad/sec. Here in Figure 6(b), the best power coefficient 
of 0.58 can be seen when the pitch angle 28 and pitch 10° is considered for the analysis for the velocity 14 
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m/s and 9 m/s. Further when the pitch angle is taken as 15° and wind speed as 11 m/s, the power coefficient 
comes as 0.52, which is near to Betz's limit and gives a promising output. 











4.2 a Picho ' —— Pitch 0 

1.1 1 Pitch 5 ‘ 4 —— Pitch 5 
= 1.0] | a eee Esti | «Pitch 10 
© 0.9 \ —— Pitch 10 g 13 +— Pitch 15 
o 0.8] \ —— Pitch 15 = AV ; 
€ 0.7 \h\ = AYI ~~ Pitch 20 
2 0.6 \\ Pitch 20 0 W —— Pitch 25 
S 0.5 W —»— Pitch 25 E AN -Pitch 28 
5 03 ANU -= Pitch 28 $ AN 

0.2 WS of) A ASR 
a 0.1 eS 5 SON Oe KOS 

Ai cna fj “ORAS 

5 10 15 20 25 77” an” 0 on 
es Sees ANUS) Wind Speed (m/s) 
(a) (b) 


Figure 6 Power coefficient curve at different pitch angle for (a) œ = 1.57 rad/sec, (b) œ = 2.22 rad/sec 


3.4. Verifications 

The numerical examination is checked by comparing it with a normal empirical equation and is 
illustrated in (3). Variation of the velocity vector is shown in Figure 7(a) and Figure 7(b), and the verification 
is shown in Table 3. It is to be noted that an additional 1 m has been considered to account for the distance 
from the root of the blade to the hub, making the blade length 44.2 m. 


Tangential velocity =r X w (3) 
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Figure 7. Tangential velocity of the blade for (a) @=1.57 rad/sec, (b) 2.22 rad/sec 


Table 4. Tangential velocity relation 


Pitch angle Velocity P Tangential velocity 
(degree) (m/s) tà ieee) Numerical Analytical Error (%) 
5 10 1.57 69.30 69.394 0.13 
5 10 2.22 97.997 98.124 0.12 


4. CONCLUSIONS 

The rotor has a rotational speed of 1.57 rad/sec produced less torque and power as compared to 
2.22 rads/sec. For both rotational speeds, torque and power show a decreasing trend for pitch angle 25° and 
28°. When rotor rotational speed is taken as 1.57 rad/sec, the optimum pitch angle obtained is 28° at 10 m/s, 
whereas when rotor rotational speed is taken as 2.22 rad/sec the optimum pitch angle obtained is 10° and 28° 
observed at wind speed 9 m/s and 14 m/s. 

In this paper, the Pitch angle plays a vital role in determining the performance parameter of the wind 
turbines. This paper gives an idea that due to variation in pitch angle, velocity, and angular velocity also 
induces some impact on aerodynamic characteristics of the wind turbine. Wind turbine performing under 
different operating conditions reflects variation in torque and power value, which can help a researcher and 
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engineer to demarcate the difference between the best pitch angle and worst and design the wind turbine 
blade accordingly. Apart from this, due to variation in pitch angle, there is uncertainty in pressure distribution 
of the blade at a specific pitch and absolute velocity, which can lead to blade deformation, and it should be 
avoided to increase the life span of the wind turbine. This variation in pressure distribution gives an idea 
about which pitch angle and which velocity should get operated to avoid hazardous like flutter and 
deformation. 
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